biseau<-function(x)
{
  colnames(x)<- c("Espcie","Ti","Tis")
  nsp<- length(levels(x$Espcie)) # Nmero de espcies
  
  if (nsp<=5)
  {
    spp <- levels(x$Espcie)
    cor <- c("blue","red","darkgreen","orange","grey")
    par(mar=c(5,5,2,2))
    graf<-plot(c(0,100),c(0,100),xlab="Proporo da espcie na descarga da viagem (%)",
               ylab="Descargas acumuladas (%)", type="n") # Grfico vazio 
    par()
    cat("Clique no grfico para adicionar as legendas\n")
    for (i in 1:nsp)
      {
        xsp<-x[x$Espcie==spp[i],]
        xsp$C<-xsp$Tis/xsp$Ti # Percentagem da espcie s na descarga da viagem i.
        xsp$j<-as.integer(xsp$C*100) # Arredondamento das percentagens (categorias j)
        agg.j<-aggregate(xsp$Tis,list(j=xsp$j),FUN=sum) # Soma do peso da sp que entrou na categoria j (=TC)
        names(agg.j)<-c("j","TC")
        dat.bis <- data.frame(0:100,rep(0,101)) # Cria um data.frame vazio
        names(dat.bis)<-c("j","TC")
        for (z in 0:100) 
          {
            ifelse(test= nrow(subset(agg.j,j==z,TC))==0, # Completa as categorias (j) que esto sem TC por 0
                  yes= dat.bis$TC[dat.bis$j==z]<-0,
                  no= dat.bis$TC[dat.bis$j==z]<-agg.j$TC[agg.j$j==z])
          }
        dat.bis$P <- cumsum(dat.bis$TC)/sum(dat.bis$TC)*100 #Contribuio relativa cumulativa 
        lines(dat.bis$P~dat.bis$j,col=cor[i])
      }
    legend(locator(1),legend=c(spp),lty=1,col=cor[c(1:nsp)], bty="n")
    return(graf)
  }
  
  
  if (nsp>5)
  {
    cat("Utilize um nmero de at 5 espcies")
  }
  
}
  